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We showtj'p that all kinds of biasing of cosmological phase transitions produce qualitatively new type of domain 
wall networks. The biased networks consist of compact, finite size, bag-like wall structures and exhibit a generic 
instability. The surface of biased networks disappears exponentially fast after a limited period of scaling. We 
argue that fluctuations of the background make the network unstable even in the case of the "symmetric on the 
average" initial distribution. We observe that the variation in parameters of the potential, like its hight, can 
influence the lifetime of the wall network, contrary to the standard beliefs. 



I. Introduction 

Topological defects provide cosmologists with a set 
of intriguing mechanisms for structure formation 
which are quite different in nature to the standard 
inflationary paradigm. Defects form at sponta- 
neous symmetry breaking phase transitions in the 
early universe, and their subsequent field ordering 
dynamics can perturb the matter and radiation 
content of the universe, leaving characteristic sig- 
nals in both the present day matter distribution 
and the microwave background. 

The nature of the defect is determined by the 
topology of the vacuum manifold following the 
phase transition. A Z2 manifold with two dis- 
connected vacua leads to domain walls, an S^ 
manifold to cosmic strings, an S^ manifold to 
monopoles, and an S^ manifold to cosmologi- 
cal texture. Much recent attention has focused 
on the gauged cosmic string, and global texture 
scenarios. However, cosmological domain walls 
have long been considered unworkable. Zel'dovich, 
Kobzarev and Okun noticed in 1975 that the en- 
ergy density of a domain wall network will eventu- 
ally come to dominate that of matter or radiation. 
More recent attempts to save domain wall scenar- 
ios, such as the the so-called "late-time" phase 
transitions, are now known to be in conflict with 
observations of the microwave background. 

The most complete study of the dynamics of 
domain walLnetworks was given by Press, Ryden 
and SpergeEl (PRS). These authors studied the 
time evolution, during the era of matter domina- 
tion, of topological domain wall kinks in a scalar 
field with a potential energy possessing two min- 



ima, both degenerate in energy. They showed con- 
clusively that such networks rapidly evolved into 
long domain walls stretching across the universe 
whose surface area, and, hence, energy density, 
persisted for a long time. This persistence, or 
scaling behavior, led both to the relatively rapid 
domination of the energy density of the universe 
by these walls and to large distortions in the cos- 
mic microwave background. Both of these results 
are incompatible with observations. 

It did appear, therefore, that domain walls 
could not have formed in the early universe. How- 
ever, these results were based on very special as- 
sumptions about the initial conditions of the do- 
main wall network. To be precise, these authors 
initialized their networks by allowing the computer 
to randomly choose, at any point on the lattice, ei- 
ther one vacuum state or the other with each vac- 
uum weighted with equal probability. However, if, 
for some reason, the two vacua were to be given 
different, or biased, probabilities, then these con- 
clusions could be dramatically altered. Further- 
more, it has become clear in recent years that non- 
equilibrium phase transitions, which can occur in 
realistic models of the early universe, generically 
lead to a biased choice of vacuum state. There 
seems to be every reason then to restudy the case 
of the walls with the modification that the domain 
wall network be initialized using biased vacuum 
probabilities. 

II. Biased Transitions 

To illustrate the way the biased phase transitions 
can actually occur we shall discuss in some de- 



tail the post-inflationary nonthermal phase tran- 
sitions in pa_ weakly coupled, out-of-equilibrium 
scalar fielcu'El. The inflationary epoch which we 
presume to precede the radiation dominated (RD) 
epoch is here the usual large-scale inflation with 
the Hubble parameter of the order of lO^'^GeVg. 
Each scalar field which lives in the Universe ud- 
ergoing very rapid, quasi-exponential expansion 
is subject to large quantum fluctuations at all 
wavelenghts, the overall magnitude of which is 
set by < (f) >^w HfSr, where St is the dura- 
tion of the inflationary epoch and Hi the Hub- 
ble parameter during inflation, in massless case, 
and < > « — I- in the case of a field of a mass 
m. More precisely, the magnitude of fluctuations 
(power spectrum) grows like -p-, so becomes very 
large for large wavelength components of the field. 
The result is, that these large wavelength compo- 
nents correlate the observables measured for the 
field over distaces as large as the blown-up infla- 
tionary horizon, while the small short-wavelength 
components decorrelate locally the measurements, 
amounting to small-scale inhomogeneities. To see 
implications this picture has for phase transitions, 
let us assume that the field has a double-well po- 
tential, the hight of the barrier between vacua be- 
ing Vo and the mass of the field - m. Then, after 
inflation, in the RD epoch, when the causal FRW 
horizon grows larger than £ ~ — , in each causal 
volume the local field starts rolling towards the 
nearest minimum of the potential. If there would 
be no fluctuations, the initial conditions for rolling 
would be the same in all horizon volumes, and the 
field would assume the same vacuum state every- 
where. However, if the characteristic scale of fluc- 
tuations. Hi, is at least comparable to the distance 
between the minima of the potential, there is a sig- 
niificant probability that in randomly choosen vol- 
umes the fluctuations have caused different min- 
ima to become local vacuum states. Now, whether 
the probability of choosing one, say (-I-), vacuum 



"It turns out that one can produce in a very simi- 
lar way nonequilibrium transitions in fields which are 
coupled strongly enough to achieve thermal equilib- 
rium - this requires inflationary stage occuring around 
the time of the transition - hence the scale of the Hub- 
ble parameter of the inflationary deSitter epoch to be 
comparable to the scale of symmetry breaking 



over the other is 1/2 or not, it is decided by the po- 
sition of the background (long-wavelengths) part 
of the field with respect to the symmetric point in 
the field space - the point eqally distant from both 
minima. To decide where the background is lo- 
cated, one should realize that the background itself 
is also a stochastic field. In fact, let us divide the 
quantum field, including all the fluctuations born 
during inflation, into 2 parts with respect to the 
scale £: the short-wavelength part T^ consisting 
of fourier components with H^^ < X < i and the 
long-wavelength part Bg with i < X < L where L 
is the blown-up inflationary horizon. With such a 
division it is clear that Te measures the difference 
or decorrelation of the field between two points at 
distance £ apart, or equvalently between "mean" 
(in the rms sense) values of the field in two boxes of 
size £ each. On the other hand, the field Bi mea- 
sures the common part of the field, which looks 
like a constant contribution if one moves over dis- 
tances £ and smaller. Of course, also Bi is in prin- 
ciple a quantum, or - at the level of quasi-classical 
description - a stochastic field, which can be char- 
acterized by means of some probability distribu- 
tion, or through moments of the distribution, like 
the mean or standard deviation. In first approx- 
imation, we can regard the Bi field as a constant 
classical background whose value is chosen accord- 
ing to a gaussian distribution with the width equal 
to the rms value of the underlying quantum field 
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Departures from this zeroth-order approximma- 
tion will be considered in the next section. It is 
clear now, that the probability that the classical 
position of the background would correspond to 
the equal probabilities of populating both vacua 
is close to zero, and that the choice of the vacuum 
in a postinfiationary transition is in general biased, 
one of the vacua becomes more likely to be choosen 
than the other. Of course, in the region interpo- 
lating between different vacua a domain wall shall 
form with a large probability. The physics of so 
formed biased networks is summarized in the fol- 
lowing section. 

III. Evolution of the Netv^rork 

The detailed procedure which has been used to 
perform numerical simulations, and a simple an- 



alytical model which can be used to understand 
main quahtative features of the result have been 
described in referencda. The focus here is to sum- 
marize the evolution of the surface energy den- 
sity of the network of domain walls in the radi- 
ation dominated epoch (the results in the mat- 
ter dominated epoch are qualitatively identical). 
Two Dimensional Simulations. In this case we 
have studiedd the evolution of A/V as a function of 
elapsed conformal time, 77, for a number of initial 
bias probabilities between p — 0.5 and 0.4. Each 
simulation ran on a 1024 x 1024 lattice until either 
no more domain walls were found in the box, or rj 
exceeded 512 (= L/2). 

For the p — 0.5 case we recover the scaling 
properties reported in PRS. Fitting the scaling 
portion (10 < rj < 100) of the curve to the power 
law 
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we find P = -0.88 ± 0.04. 

Moving away from the p = 0.5 case, one sees 
a dramatic departure from self-similar scaling. In 
each case there is an exponential cut-off in the 
ratio A/V at some characteristic time. For the 
cases of p close to 1/2, that is for p=0.457 and 
0.45, we find that the curves are well fitted by a 
function of the form 
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However, for the cases of p near 1 — pc = 0.407, 
that is, p — 0.425 and 0.4, a simple exponential 
suffices: 
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Values for D and fj, averaged over 5 runs for each 
value of p, are given in Table |j 

For p close to 0.5 the domain wall network ap- 
pears to enter a quasi-scaling regime in which A/V 
scales ex ry"^ , before eventually being exponentially 
cut-off a.t T] = fj. As p ^ 1 — Pc = 0.407, how- 
ever, fj rapidly approaches the resolution size of 
the grid, and no evidence of early scaling is seen. 
This behavior continues as p drops below the criti- 
cal threshold, as indicated in the p = 0.4 example. 
We find that for the cases oi p < 0.4 the exponen- 
tial cut-off of A/V becomes even more precipitous, 
with the characteristic scale fj rapidly approaching 
the lattice grid size. 



To conclude, in the two-dimensional simula- 
tions we see persistent scaling behavior precisely 
a.t p = 0.5. For p below 0.5 but above the critical 
threshold l—pc = 0.407, we see scaling for a finite 
time which is then exponentially cut-off at some 
conformal time fj. The value of fy, which becomes 
very large as p -^ 0.5, decreases rapidly as p ap- 
proaches the critical threshold. Near p = 0.407, 
and below it, no scaling behavior is seen and the 
behavior is well described by a simple exponential 
for all conformal time. 
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TABLE I.; Fits to the plots of AjV against r] for different 

initial bias probabilities,^ p, in too dimensions, using the 

functional forms (H) and (W) given in the text. 

Three Dimensional Simulations. The three di- 
mensional simulations are run on a 128'^ grid. 
Again, the self-similar evolution seen in PRS for 
the p = 0.5 case is reproduced well in the time 
range 2wo < rj < L/2. Measuring the logarithmic 
slope of A/V versus rj between these times we find 
u= -0.89 ±0.06. 

With p ^ 0.5, we see qualitative features in 
the three dimensional runs similar to those of the 
two dimensional cases. The one major difference 
is that the turn-over into an exponential decay is 
seen to occur much earlier, so much so in fact that 
all the cases considered here {p = 0.47, 0.48, 0.49) 
are well fitted by a simple exponential curve with 
no initial pseudo-scaling regime. Using a fitting 
function of the form of equation flthe values of 
the fj for each p are given in Table ||. 
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TABLE IL: Fits to the plots of A/V against rj for different 

initial bias probabilities, p, in three dimensions using the 

functional form (kj) given in the text. 



To conclude, in the three dhiiensional case we 
see persistent scaUng behavior precisely at p = 0.5. 
For any value of p below 0.5, little or no scal- 
ing behavior is seen and the behavior is well de- 
scribed by a simple exponentiaLfor all conformal 
time. Background Fluctuationn. The new re- 
sult which we want to discuss here is the influ- 
ence of the fluctuations of the background on the 
behaviour of the network. In the simulations de- 
scribed so far the background value of the stochas- 
tic field was assumed to be a constant randomly 
located in the field space (the freedom of choosing 
this background value is reflected in the freedom of 
choosing the bias parameter p). However, we know 
from the previous section that in fact the back- 
ground is not strictly homogeneous, as there must 
be fluctuations around the zeroth-order approxim- 
mation constant background value. In the simu- 
lations whose results are depicted in Figure we 
are simulating the fluctuations of the background 
through fluctuations of the probability p. At each 
lattice site we draw p according to the gaussian 
distribution with some central value p — Po and 
the width equal to cr = 0.25po- As one can see 
from the plot, this has a dramatic impact on the 
bahaviour of the equilibrium network with previ- 
ously constant value p = Po = 0.5. Now, with 
fluctuating p, the network becomes unstable, and 
its surface again decays quasi exponentialljcl. In- 
tuitively one can understand this result noticing, 
that the average value of the probability over a 
randomly choosen small subvolume of the lattice 
is typically different from po — 0.5, hence each 
small part of the network behaves like a biased 
network, and as these subnetworks are unstable, 
the averaging over many subnetworks cannot pro- 
duce the equilibrium bahaviour, despite the fact 
that the sub-probabilities average to po = 0.5. 

The other case shown in the picture (dotted 
curves, p = 0.3) illustrates the general tendency 
of biased networks to persist for a longer period 
of time (slower decay rate) when the probability 
(in fact - the background field configuration) fiuc- 
tuates between different lattice sites. The fluctu- 
ating p — 0.5 case shows that in reality all the 
networks formed during actual cosmic transitions 
will be biased, hence after a finite scaling period, 
typically rather short, will disappear dissipating 
their surface energy in the form of scalar (maybe 
also gravitational) waves. 
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FIG. 1.: The evolution of the surface area of the network 
when background fluctuations are present. The thin solid 
line corresponds to the constant value of p = 0.5 while 
the thick solid line shows the situation where p fluctuates 
around po = 0.5. Instability in the later case is clearly 
visible. The thin dotted line presents the evolution of the 
network created with constant value of p = 0.3. The thick 
dotted line gives evolution of the biased network when p 
fluctuates around po = 0.3. All curves correspond to d=2 
dimensions. The verical axis gives surface of the network, 
the horizontal axis gives conformal time. 



Wall ThicknesSi'a. The next observation which 
can have some significance for the problem of the 
possible impact of biased networks on structure 
formation is that the values of 9 and fj were found 
to be weakly dependent on the value of wq, the do- 
main wall thickness. Changing the wall thickness 
affects the network evolution through two closely 
balanced effects. Increasing the height of the po- 
tential barrier between the two vacua, Vo — 5^, 
makes transition from one vacuum to the other 
energetically less favorable. But, this is compen- 
sated somewhat, not necessarily exactly, by the 

1 ^ 
increased domain wall surface tension, cr = t^-^— . 

2 Wo 

However, for a quartic potential, the value of Vq 
clearly does not scale out of the equations of mo- 
tion, so this balance need not be exact. (In the 
domain wall scenario, the scalar field is off the vac- 
uum manifold at many points throughout space, 
making the evolution more dependent on the rel- 
ative values of Vb and ipQ than, say, a texture sce- 
nario, where the non-linear sigma model works 
well, and evolution is effectively independent of 
the potential barrier between vacua.) This is il- 
lustrated in the Figure |[ 




FIG. 2.: The curves show the dependence of the evolution 
of the network on the parameters of the potential through 
the domain wall thickness Wo- There is a tendency for the 
network to persist for a longer period of time for thinner 
walls. The figure shows the case of p = 4 in d=2 dimen- 
sions. The results in d=3 are qualitatively the same. The 
vcrical axis gives surface of the network, the horizontal axis 
gives conformal time. 



IV. Conclusions 

From the presented results it becomes clear that 
biased transitions solve the domain walls problem 
very efficiently and naturally due to the generic in- 
stability of the network with respect to the depar- 
ture from p = 0.5 and with respect to fluctuations 
of the background. 

The interesting and important for structure 
formation problem which remains to be clarified 
is whether there is a possibility of having unsta- 
ble networks with scaling periods long enough to 
show up on the — curve. This would be an inter- 
esting possibility in the context of hot dark matter 
scenarios. The other point which has to be consid- 
ered more detaily is interaction of biased networks 
with matter. 

Recenly we became aware of the work of S. 
Larsson, S. Sarkar and P. Whitra which confirms 
and extends our invstigations of biased networks. 
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